Elastic and non-linear stiffness of graphene: a simple approach 
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The recent experiment [Science 321, 385 (2008)] on the Young's modulus and third-order elastic 
stiffness of graphene are well explained in a very simple approach, where the graphene is described 
by a simplified system and the force constant for the non-linear interaction is estimated from the 
Tersoff-Brenner potential. 
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Graphene has received a large number of interests since 
its discoveryji*^ Among others, a recent experiment has 
shown that the graphene has excellent mechnical proper- 
ties with very large Young's modulus;^ This experiment 
also found obvious non-linear effect for the graphene in 
the large strain regime. Theoretically, the elastic prop- 
erties in the graphene can be studied in the continuum 
mechanics approach^ or the ab initio method.^ In Ref. 0, 
Pugno derived the third-order elastic stiffness (TOES) by 
relating it to the coefficient of thermal expansion- . In this 
paper, we describe the graphene in a simplified system 
with very simple interaction potential. There is no vari- 
able parameters in the potential, where the non-linear in- 
teraction can be deduced from the Tersoff-Brenner (TB) 
potentiali^i^ The recent experiment for the Young's mod- 
ulus and the TOES of the graphene are well explained by 
our model. 

Graphene samples in the mechanical experiment^ usu- 
ally have a radius larger than 0.75 /im. In this large 
radius two dimensional sheet, the number of atoms on 
the boundary (TVb) is much smaller than that of the in- 
ner atoms (Ni). It can be estimated as: 
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where r is the radius of sample, and b = 1.42 A is the 
C-C bond length in the graphene. s — -^(\/36)^ is the 
area of the unit cell in the graphene. The ratio of these 
two numbers is: 

= 1.5V36/r w 5 X 10-^ 

where r — 0.75/im is used. It shows that the number 
of boundary atoms in the sample is about four orders 
lesser than the inner atoms. As a result, the contribution 
of the boundary atoms to the total energy is also four 
orders smaller than the inner atoms. In this sense, we 
can ignore the contribution of the boundary atoms to 
the total energy in the graphene. In the following we 
consider only the inner atoms. 

Because of the translational symmetry, all unit cells 
in the graphene are equivalent to each other. We can 



FIG. 1: (Color onhne) Configuration for the simphfied sys- 
tem, (a) is the equilibrium position, (b) is the configuration 
with strain. The inner atoms 4 and 5 are optimized into new 
equihbrium position 4' and 5' with four outside atoms fixed. 



consider only one unit cell as a representative of the 
graphene, as shown in Fig. [Ha), where atoms 4 and 5 
are the two non-equivalent carbon atoms in the repre- 
sentative unit cell. 

In this paper, the direction with angle 9 — Q directs 
from atom 4 to 5 in Fig.[IJa). And 6* = 7r/2 is the vertical 
direction which is named armchair direction throughout 
this paper. 

The interaction potential we used includes both lin- 
ear and non- linear terms. The linear interaction is the 
frequently used bond stretching interaction V/: 
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where, b (bo) is the strained (unstrained) C-C bond 
length in graphene. The force constant ki — 305Nm~-'^ 
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TABLE I: Parameters in the Tersoff-Brenner potential. See text for the meaning of each parameters. 
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is taken from Ref. llCt where this potential was appHed 
successfully to explain phonon properties in the carbon 
nanotubes and graphene layers. 

The non-linear interaction we applied has the form: 



(2) 



and the constant kni can be evaluated from the commonly 
used TB potential. Similar expansion has also been done 
in Refs. Inllil. 

In the TB potential, the energy is expressed as: 

Vs(r) = VRir)-B,,-VAir). 

Vr and Va are the repulsive and attractive energy: 
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with the cut-off function /c(r 
fcir) = {\[ 
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The many-body coupling parameter is: 
1 
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The angle function G{6ijk) is: 



G{9ijk) — oo 



df, dl + {1 + cosie,jk)Y 



All parameters in the TB potential are listed in Table. 
The ratio of kni / ki can be estimated from the TB poten- 
tial as following. 

(1). In the equilibrium structure of the graphene, all 
of the angles equal to |7r. So the angle function G{6ijk) 
can be simplified: 



G{9) = ao 



1 + ^ 



dl d2 + (l + cos(|7r))' 




FIG. 2: (Color online) The strain energy v.s strain. Strain 
energy at e = —5% is about 15% higher than that at e = 5% 
implicating the non-linear interaction. The fitted function is 
E{e) = iFe^ + i_De^ with Y = 1.131 TPa and D = -2.360 
TPa. 



(2). Due to the small value of G{9), we have: 
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(3) . Because ii^^^=2 A in the cut-off function is about 
40% larger than the equilibrium bond length 1.42 A in 
the graphene, we simply set: 

fcir) « 1. 

(4) . We expand the exponent function in the repulsive 
and attractive energy in terms of r — i?'*^-* : 
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with X = V2S(3ir - R'^"^). 

(5). Finally, the total binding energy is: 

Vsir) « VRir)~VAir) 
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FIG. 3: Young's modulus of the graphene. 



FIG. 4: The thhd-order stiffness D in the graphene. 
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where Vq = Z^'*^^ . The vanish of the 1st order term in this 
final expression indicates that our above approximations 
are physically correct. 

(6). Applying parameters in Table. (P), we obtain the 
value of kni/kf. 

km 1 S'+l r— 

So, the non-linear constant is /c„; = 915Nm~-'^A~-'^. 

In the mechanical experiments on the graphene, the 
sample is stretched by force and the boundary is then 
fixed to measure the relation between the force and 
the strain. Similarly, as shown in Fig. [TJ the sys- 
tem is strained (with strain e) in a particular direction 
(Fig. [11(b)). Then the outer four atoms 0, 1, 2, and 3 are 
fixed, which simulates the experimental conditions. The 
two inner atoms 4 and 5 can move freely to achieve a 
new configuration 4' and 5' with minimum energy E[e) . 
This optimized configuration is the equilibrium position 
for the system under strain. 

Fig. [5] shows the relation between the density of strain 
energy and the strain. The strain in this figure is added 
along 6 = 7r/2 direction with the value in [— 5%,5%], 
which is also the magnitude of strain added in the 
experimenti^ The non-linear effect can be clearly seen 
from this figure, since £^(—5%) is about 15% higher 
than E{5%). Fitting this curve by function E{e) = 
^Ye^ + ^De^ , we can get the Young's modulus Y = 1.131 
TPa and the TOES D = -2.360 TPa. 

Fig. [3] is the dependence of the Young's modulus on 
the direction angle 9. The average value over 6 is about 
0.83 TPa. We assume that the experimental measured 



Young's modulus value 1.0 ± 0.1 TPa is an average over 
9. This assumption can be qualitatively understood in 
terms that the atomic force microscope tip touches the 
center of the circular graphene sample; thus generating 
strains equally in all radial directions. Our calculation 
is in good agreement with this experimental result. As 
shown in Fig. IH the average value for TOES D is about 
-1.3 TPa, which is comparable with the corresponding 
experimental value -2 ±0.4 TPa. Our result is a lit- 
tle smaller due to the simplicity of our model, which 
may underestimate the non-linear interaction in the sys- 
tem. The main reason for the underestimation is that 
the interaction V/ can dominate the total interaction in 
graphene. However, there are some other weaker inter- 
actions in the graphene, such as twisting interaction. As 
these weaker interactions are neglected in present model, 
the Vi with force constant ki underestimates the linear 
interaction. Consequently, the non-linear force constant 
from kni = 3ki will underestimate the non-linear interac- 
tion in the system. 

In conclusion, we calculate the Young's modulus and 
TOES of the graphene with a simplified system and very 
simple interaction without any variable parameters and 
explain the recent experimental results nicely. 

We further remark that a large piece of graphene can 
be regarded as a thin plate. In a plate, the non-linear 
effect of the third order arises from terms in the elas- 
tic energy which are cubic in the strains This non- 
linear effect leads to a non-linear equation of motion for 
the system, which indicates the coupling between dif- 
ferent phonon modes. In the graphene, the non-linear 
interaction has exhibited itself in different phenomena. 
As a direct result of the non-linear interaction, the Ra- 
man G mode shows a red-shift with the increase of 
temperaturcj^ Very recently, the non-linear interaction 
in carbon nanotubes is confirmed to be the origin of the 
intrinsic localized mode, which leads to the Stone- Wales 
defect under axial tension. This effect should exist in 
the graphene. 
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